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I OBJECTIVE (U) 


(S/CL-3/NF) The objective of this effort is to investigate the 


possible effects of ambient geophysical/ low-frequency electromagnetic 


factors on remote viewing (RV)* performance as a potential aid to 


increasing the performance levels of Army INSCOM remote viewers. 


* 
(U) RV (remote viewing) is the acquisition and description, by mental 
means, of information blocked from ordinary perception by distance 
or shielding. 
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II INTRODUCTION (U) 


(S/CL-3/NF) SRI International is tasked to conduct a study for 
Army INSCOM to investigate a potential correlation between remote viewing 
(RV) performance and ambient geophysical/extremely-low-frequency electro- 
magnetic (ELF) activity. The possibility of such correlation is indicated, 


for example, by studies showing psychophysiological effects*’?" 


3,4 


and 
behavioral changes associated with ELF electro-magnetic fields. The 
geophysical variables of interest include such factors as ELF intensity/ 
fluctuations, ionospheric conditions, geomagnetic indices, sunspot number, 
and solar emissions (e.g., X rays and solar flares). The questions to be 
answered in this program are 

¢ Do geophysical/performance correlations exist such that 


measurement of the ambient geophysical variables could 
be used as an indicator of expected performance? 


e If so, can optimum performance windows be identified? 


(U) The structure of the program that will address the above issues 


consists of 


e A literature search. 
e Real-time ELF measurements 


— SRI (Menlo Park, California location) 
— Time Research Institute (Los Altos, California 
field station). 


e Real-time geophysical data acquisition via NOAA 
(National Oceanic and Atmospheric Administration) 
Westar IV satellite downlink. 


e Computer correlation studies of RV performance versus 
variables of interest. 


* 
(U) References are listed at the end of this report. 
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III PRESENT STATUS (U) 


A. (U) General 


(U) In order to accomplish the goals set out in Section II, the 
program has been designed as a joint effort between SRI International and 
Time Research Institute of Los Altos, California, with SRI as the prime 
contractor. Time Research Institute is a research organization that 
specializes in temporal analysis of geophysical variables and their po- 
tential correlation with phenomena of interest, such as weather patterns, 


earthquakes, and so forth. 


(U) With regard to the present effort, Time Research Institute is 
responsible for establishing the appropriate hardware and software systems 
for collecting and analyzing data concerning environmental conditions and 
their correlation with RV performance. The purpose of the correlation 
study is to determine whether RV performance is enhanced or degraded by 
measurable changes occurring in the geophysical (including solar- 
terrestrial) environments. The specific data bases under consideration 


in this effort are given in Table 1. 


(U) Should correlations between geophysical variables and RV 


performance be found, the application potential of the effort is twofold: 


(1) Time periods in which enhanced RV performance might be 
expected could be identified, resulting in increased 
quality and accuracy of information obtained through 
such channels; similarly, time periods in which degraded 
RV performance might be expected could be avoided. Thus, 
optimum performance windows would be identified. 


(2) An increased understanding of the types of environmental 
changes that correlate with RV performance could provide 
clues as to the mechanisms involved in RV functioning. 
Such knowledge would lead to more focussed research on 
factors that could enhance RV performance, and would 
also provide information critical to the development of 
defensive countermeasures against RV. 


5 
Approved For Release Pag RB Ass STF TE TI9®00788R001 800160001-6 


Approved For Release 2009 REEL: A'S SPPIE TG /88RO0 18001600016 


Table 1 


(U) GEOPHYSICAL DATA BASES 


Solar-Terrestrial 
— Geomagnetic 


Ground-measured indices Ap, K, etc. 
Satellite-measured intensities 


Solar flux (MHz) 


15,400 1,415 
8,800 606 
4,995 410 
2,800 245 
2,695 


Sunspot number 

Solar flares 

Interplanetary magnetic field 
Solar wind (Pioneer XII) 
Protons 


Cosmic ray indices (neutron monitor) 


Ionospheric Measurements 


— Sudden ionospheric disturbances (SIDS) 


— Auroral electrojet 


— Radio propagation quality indices 


ULF/ELF 
— 30 frequencies (from 1 to 30 Hz) 


UNCLASSIFIED 


(U) The tasks Listed in Section II (literature search, real-time 
ELF measurements, real-time geophysical data acquisition via satellite 
downlink, and correlation studies of RV performance versus geophysical 
variables of interest) have been prioritized with the goal of producing 
the longest possible ELF data base during the period of this contract. 
Therefore, while all of the tasks are being pursued in parallel, ELF- 


related tasks have been the focus of attention to date. 
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B. (U) Literature Search 


(U) The purpose of the literature search is to acquire and 
integrate information from recent work done in the field of biological 
response to ELF and geophysical conditions. Approximately 30 keywords 
are in use as input to a computerized literature search. Forty manu- 
scripts have been obtained to date and are in the process of being 
reviewed, plus additional sources of literature have been identified and 


will be retrieved as priorities permit. 


Cc. (U) ELF Measurements 
1. (U) Introduction 


(U) Although the ELF frequency range (3 to 300 Hz) has been 
studied in some detail, many unknowns remain. For example, although it 
is known that ELF frequencies generated by geophysical means (e.g., elec- 
trical storm activity) tend to distribute themselves globally, little 
information is available on the variation of the ELF environment from 
location to location. Therefore, local variations may exist that are 
caused by both manmade sources, and by the geological structure of the 
area. In the San Francisco Bay Area, manmade sources that generate ELF 
on a local scale include motors, telephone lines, power lines, and 
electrical subways [Bay Area Rapid Transit (BART)], and it needs to be 
determined whether the emission from such sources constitutes a signif- 


icant contribution to the omnipresent global ELF field. 


(U) In order to address the above issue, two ELF monitoring 
stations are being set up--one at SRI Menlo Park (in the RV Laboratory) , 
the other at the Time Research Institute field station, 17 km distant. 
It is anticipated that the SRI environment may be a "noisy" one due to 
the large amount of electronics known to be in the area. Data from the 
two sites, taking the field site as a reference, are to be compared in 
order to begin to differentiate the naturally-occurring ELF from the 


manmade noise occurring at the location where RV is being carried out. 
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2. (U) Los Altos Site 


(U) Since May 1982, Time Research Institute has been operating 
an ELF monitoring site in Los Altos, collecting data twice daily for the 
purpose of correlating ELF disturbances with various phenomena of interest. 
In this period, analysis techniques were developed that are directly 


applicable to the present task. 


(U) One of the first tasks was the upgrading of the Los Altos 
ELF monitoring site to provide coverage during power interrupts. Research 
was performed to determine the best power-interrupt system, including 
generators. The selection criteria chosen for the generator sought to 
optimize power output, cost effectiveness, and papubation for reliability, 
with the consideration that this system might serve as a model for addi- 
tional sites in the client community. On this basis, a 3500-W Kubota 
generator was bought and installed at the site. A PTI Industries "Data- 
shield" device was also purchased and installed, for use in conjunction 
with the generator (which must be started manually). This device powers 
the required electronics for a period of twenty minutes on its own while 
awaiting generator startup. Furthermore, an automatic alarm telephone 
dialing system dials up as many as four individuals should a power failure 
occur while no one is at the site. The two devices working together have 
protected the system on numerous occasions from power interruptions-- 
including an 8-hr outage planned by Pacific Gas and Electric Company. 
Thus, since installation of the power-interrupt equipment, there has been 


no loss of data collection/storage. 


3. (U) ELF Data Acquisition Systems 
a. (U) Basic System Design 


(U) With the requirement that two ELF monitoring sites 
be implemented for the program--one at SRI and one at Time Research 
Institute--it was decided that the two systems would be made identical. 
In this way, differences between the two systems would be minimized, thus 
reducing the opportunity for artifactual differences between the two sys- 


tem outputs. 
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(U) Figure 1 is a block diagram of the basic ELF data 
acquisition system. The ELF signal is collected by an antenna, amplified, 
and then digitized by an analog-to-digital (A/D) converter so that the 
signal can be input into an IBM PC Jr. computer for the purpose of analy- 
sis by a fast-fourier-transform (FFT) program on at least an hourly basis. 
The recorded data are then transferred by floppy diskette to an IBM XT 


computer for further handling. 


AMPLIFIER 
A/D 
CONVERTER 
COMPUTER 
IBM PC JR. 
Re ee 


HOURLY 
DATA 


UNCLASSIFIED 


INCREASES SIGNAL STRENGTH © 


MAKES SIGNAL READABLE BY 
COMPUTER - 


DETERMINES INTENSITY OF SIGNAL 
FREQUENCY COMPONENTS 


FLOPPY DISKETTES TRANSPORTED 
TO IBM XT COMPUTER 


FIGURE 1 (U) ELF DATA-ACQUISITION SYSTEM 


(U) As indicated in the above system description, an 
integral part of the data acquisition system is computerized record 
keeping, using IBM systems--both the IBM PC Jr. and the IBM XT. The 
software is written and compiled on the XT, used as a master system, and 
then run on the PC Jr. (the PC Jr. is not itself large enough to run a 
compiler, nor are there compilers written for it). Beyond this, however, 


in spite of the much publicized "compatibility" between the various 
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IBM PC systems, some development time had to be expended to ensure that 
programs compiled on the XT (that were concerned with communication with 
external devices) could be run on the PC Jr.--because they handle the 


addressing of communication ports differently. 


b. (U) FFT Program 


(U) A consultant was hired from the Radioscience Labora- 
tory at Stanford University to critique the ELF system at Time Research 
Institute as it stood at program start. As part of that critique, he 
recommended changes in the FFT software to increase its accuracy and 
running speed. As a result, new FFT software was written in compiled 
BASIC and then debugged. The program description is given in Appendix A. 
This task was completed in mid-April. 


C. (U) System Electronics 


(U) The prototype system of Figure 1 was assembled and 
installed for our purposes at the Los Altos site. Early tests indicated 
that modification of the existing electronics was required because 
(1) the system was sensitive to ground-Loop problems, thus the preamplifier 
had to be redesigned to include an isolation amplifier, and (2) when the 
new system was installed in mid-June, it was found that a slowly-varying 
dc level was superimposed on the incoming ELF signal, resulting in ex- 
cursions that exceeded the limits of the A/D converter at amplification 
levels required for good signal analysis. Therefore, new circuitry was 
designed to eliminate the dc problem. With these changes, the ELF detec- 
tion system is scheduled for completion of testing, debugging, and 
calibration at the Los Altos site in August, before its sister system 
is installed at SRI. 


d. (U) ELF Electronics/Software Subtasks (Status) 


(U) The status of the ELF electronics/software subtasks 


is as follows: 
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(U) 


Subject to the requirement that the basic microprocessor 
units to be used in the program would (1) be able to com- 
municate with an IBM XT, (2) have at least 64K of memory, 
(3) possess diskette storage capability, and (4) be cost 
effective, the IBM PC Jr. was selected from among the 
various alternatives, and two units were purchased. 


Preamplifiers, low bandpass filters, and amplifiers have 
been designed, and one complete prototype system has been 
assembled. Amplifiers and filters have been assembled 
for the second system to be installed at SRI. 


Design specifications have been completed for the ELF 
antenna. Assembly of this antenna is pending the results 
of a calibration task (described in Section III.D.3.e). 


A survey of A/D converters that would be compatible with 
the other system components was completed, and the 
selected units were purchased. 


Communication between the IBM PC Jr. and the A/D con- 
verter has been established, enabling the computer to 
read the incoming signal. 


Software has been developed and debugged that: (1) reads 
the communications RS 232 port of the PC Jr. input from 

the A/D converter, (2) performs FFT analysis of the signal, 
and then (3) writes half-hourly averages of 19 different 
frequencies (from 1 to 29 Hz) to a computer diskette. 
Further software refinements will continue to be made, such 
as determining daily maximum values for each frequency. 
This software has been implemented and debugged. The sys- 
tem is ready to begin data acquisition upon implementation 
of the dc-level-controlling hardware. Minor enhancements 
of the software will continue. 


The Los Altos system is now in operation and is being 
tested. It is anticipated that the SRI system will be 
implemented in August. 


e. (U) ELF System Calibration 


(U) System calibration is proceeding. The Stanford 


University consultant mentioned earlier is an expert in the areas of ELF 


and VLF measurement, antenna design, and spectral analysis. His cali- 


bration program is being carried out at the present Los Altos site using 


specialized, sophisticated instruments from the Stanford Radioscience 


Laboratory. 
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(U) As part of the calibration process, certain power- 
line and power-supply "noise" sources were identified that could have 
produced artifacts in the historical and ongoing data. | Therefore, in- 
creased electronic-filtering measures are being incorporated into the 
system, including software changes to recalibrate past data already in 


storage. 


(U) Finally, as part of the calibration process, three 
systems are to be monitored simultaneously: (1) the original unmodified 
Los Altos system on which the historical ELF data have been acquired; 

(2) a new system (described in this report) using an ELF wire-coil 
antenna designed by and on loan from the Radioscience Laboratory at 
Stanford University; (3) the new system, but with a higaneenia (oak tree) 
as an antenna, which is a procedure recommended by Stanford (see Nature 
reprint, Appendix B). The results from the calibration tests will be 
used to calibrate the historical data base, and to fix the final antenna 


design. 


D. (U) Satellite Downlink Geophysical Data-Acquisititon System 


(U) A near real-time satellite downlink system for solar-terrestrial 
data has recently become available from NOAA (National Oceanic and Atmo- 
spheric Administration). With this unit, it is possible to provide 
immediate feedback and/or analysis in conjunction with RV sessions. 
(Normally, there are long delays in procuring solar-terrestrial data; 
without the downlink, delays of 10 days to 6 months are standard.) The 
downlink system provides for accumulation of a detailed data base directly 


on computer diskettes. (See Appendix C for an item-by-item description.) 


(U) A satellite controller and a dish antenna for the downlink 
system were ordered and installed at the Los Altos site early in the 
project (see Figure 2). At the time of this writing, specifications for 
data-acquisition software for the IBM PC Jr. have been completed, and 
first-order software has been written that captures the data to computer 
diskettes. Because of the large volume of data transmitted each minute, 


a double-density, double-sided diskette fills in about 21 hours. Disk- 
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ANTENNA 


CONTROLLER 


IBM PC JR. 


<a 
DATA 
DISKETTE 


IBM XT 
COMPUTER 


BURROUGHS B6900 
COMPUTER 


UNCLASSIFIED 


FIGURE 2 (U) REAL-TIME GEOPHYSICAL DATA 
ACQUISITION VIA WESTAR IV 
DOWNLINK 


(U) 


ettes have been changed on an almost daily basis since early April in 


order to begin to acquire a long-term data base. 


(U) To increase the amount of data that can be stored on a 
diskette, it is necessary to separate numerical data of interest. The 
necessary software is now being written that will be able to identify 
the various data types as they are transmitted, so that only the data of 
interest will be transferred to an appropriate data file on the diskette. 


Statistical analysis can then be done on the data in this form. 


(U) Other data bases are continuing to be maintained for this 
project by Time Research Institute. Files of 2800-MHz solar flux, the 
planetary magnetic activity index (Ap), the Anchorage magnetic index, 
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(U) 


and the Stanford mean solar magnetic field are updated on a weekly or 
monthly basis as the data are available. ELF data from the "old" system 
are recorded twice daily; ELF data from the "new" system are being 

recorded at the following intervals: half-hourly averages (48 times a day), 
and two sets of half-daily averages (twice a day) at 00:00 and 12:00 UT, 

and at 04:30 and 16:30 UT. Whole-day averages are also being recorded. 


(U) To summarize the status of the satellite-downlink data- 
acquisition system: 
¢ The downlink-geophysical system is in place and in operation 
at the Los Altos site. 
e Partial data acquisition from the downlink is in place. 


e Software for the final data-acquisition system will be 
completed and implemented in August. 


E. (U) Geophysical Data/RV Performance Correlation Analysis 


(U) Data for RV sessions are to be analyzed for statistical cor- 
relation with respect to the battery of geophysical data sets listed in 
Table 1, and those acquired by Time Research Institute via measurement 
(Figure 3). Recording of RV and geophysical data is now in progress. 

The overall system for data acquisition and analysis of RV performance/ 
geophysical data is depicted in Figure 4. When enough data have been 
collected toward the end of the contract period, analysis will be performed. 
The tasks described in earlier sections are in preparation for this task, 
and therefore have received the bulk of the effort. Certain subtasks in 
the analysis task, which require a longer leadtime, have, however, already 


been completed in preparation for the analysis. 


(U) The primary statistical program that will be used to scan the 
data for possible relationships is called EPOCH ANALYSIS. This program 
reads two files simultaneously. The first file is an event file, the 
second a data file. The program first reads an event, then scans the 
data temporally backward and forward in time around the event. This 


information is stored, a second event is read, and so forth. When all 
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LOS ALTOS STATION SRI STATION 


NTT, 
WESTAR IV DOWNLINK ELF DATA 
IBM PC JR. No. 1 IBM PC JR. No. 2 
ane 


——e- 
DATA 


DATA 
DISKETTE 


DISKETTE 
IBM XT COMPUTER 
BURROUGHS 86900 COMPUTER 


FIGURE 3. (U) REAL-TIME GEOPHYSICAL DATA-ACQUISITION SYSTEM 


ELF DATA 


IBM PC JR. No. 3 


DATA 
DISKETTE 


UNCLASSIFIED 


(U) 


the events and surrounding data have been read, a printout is created 
that lists appropriate cross-correlation statistics between event and 


data elements. 


(U) Preliminary scans of data generated during an approximate 
100-site series with one remote viewer have been carried out. The 
session quality was graded on a scale of 0 to a and correlations 
between solar magnetic field and solar sunspot number were investigated. 
Some correlation between RV performance and solar sunspot number was 
found, which, if substantiated by further data, would indicate the 
possibility that performance might improve immediately after a peak in 
the sunspot number, and would deteriorate just before the sunspot number 
peaks in its 27-to-29-day cycle. This result is based on data points 
that are too small in number to be taken seriously at the point, however, 
and is mentioned only to give an example of the types of correlations 


that will be sought out and examined during the course of study. 
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(GEOPHYSICAL DATA) 


PERFORN- 
ANCE 


PHYSICAL 
DATA 


IBM XT COMPUTER 
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ANALYSIS 
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FIGURE 4 (U) GEOPHYSICAL/PERFORMANCE DATA-ANALYSIS SYSTEM 
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IV SUMMARY (UV) 


(U) Approximately 80 percent of the project's data-acquisition 
tasks have been completed, and all of the equipment and hardware have 


been purchased and delivered. 


(U) System calibration should be completed in the near future, 
and ELF and downlink-data acquisition will have begun in their final 
formats. 

(U) At the above point, the focus of effort will turn to analysis 


of past and present geophysical and ELF data, soon to be followed by the 


initiation of correlation studies of these data against RV performance. 
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Appendix A 


FAST FOURIER TRANSFORM ROUTINE FOR ELF DATA (U) 
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FAST FOURIER ROUTINE FOR ELF DATA (U) 
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 ehahaheket ahahaha heheheh eRahelelshahehehs eek ehehehehehedahehehahehelshebehehehabehehshehehelehehelielshohshehehihshehehelehshohehshehelehehal 


Pm * 
iad PROGRAM DESCRIPTION w 
1" * 


gH ee Be ee Ww He eee ee ee We ee WO ee ee 


i] 

sTHIS PROGRAM USES THE FAST FOURIER TRANSFORM (FFT) ALGORITHM TO 
yCALCULATE THE SPECTRAL MAGNITUDES OF AN ARRAY OF CLOSE-PACKED REAL INPUT 
;DATA POINTS. THE PROGRAM IS ORGANIZED AS FOLLOWS: 

3 

1. AT THE START, DATA POINTS ARE STORED IN THE REAL ARRAY X(14), 
WHERE + RUNS FROM @ TO NX-1. NX MUST BE AN INTEGRAL POWER OF 2, 
OF THE FORM NX = 24MX (IE, MX = LOG2(NX)). SUCCESSIVE ELEMENTS 
OF X REPRESENT SUCCESSIVE SAMPLES OF AN INPUT SIGNAL, SAMPLED AT 
REGULAR INTERVALS OF TIME DT. 


2. THE REAL ARRAY XC i) IS TREATED FOR THE FFT AS A COMPLEX ARRAY 
OF NX/2 PAIRS OF REAL AND IMAGINARY ELEMENTS. THAT IS, THE REAL 
ELEMENTS OF THE ARRAY ARE THE EVEN-NUMBERED INPUT SAMPLES, ANDO 
THE IMAGINARY ELEMENTS OF THE ARRAY ARE THE ODD-NUMBERED INPUT 
SAMPLES. AFTER THE FFT IS CALCULATED, AN ADDITIONAL STEP IS USED 
TO EXTRACT THE SPECTRUM OF THE REAL INPUT DATA. TRANSFORMING 
CLOSE-PACKED DATA IN THIS WAY, EVEN THOUGH IT REQUIRES AN 
ADDITIONAL STEP FOR THE REAL TRANSFORM EXTRACTION, IS FASTER THAN 
TRANSFORMING A COMPLEX ARRAY OF NX REAL-IMAGINARY PAIRS (2*NX 
ELEMENTS), WHERE ALL OF THE INITIAL IMAGINARY VALUES ARE ZERO. 

t 


3. THE FFT 1S CALCULATED AS FOLLOWS: 


A. THE INPUT PAIRS IN X ARE REARRANGED IN BIT-REVERSED 
ORDER. 


B. THE FFT IS CALCULATED USING DECIMATION IN TIME, WITH 
ANGLE ARGUMENTS IN EACH SUB-DFT APPEARING IN NATURAL 
(IE, INCREASING) ORDER. 


C. AFTER THE LAST PASS THROUGH THE FFT ALGORITHM, THE 
ARRAY X CONTAINS SPECTRAL VALUES IN NORMAL ORDER, WITH 
EACH EVEN POINT A REAL VALUE AND EACH SUBSEQUENT ODD 
POINT THE CORRESPONDING IMAGINARY VALUE. 


4. FOLLOWING THE FFT THE REAL TRANSFORM IS EXTRACTED. NOTE THAT 
IF THE SAMPLING INTERVAL (THE TIME BETWEEN SUCCESSIVE SAMPLES) IS 
DT, THEN THE TOTAL SIGNAL INTERVAL PROCESSED IS NX*DT (THE FFT 
ASSUMES THAT THE MISSING RIGHT END POINT IS THE SAME AS THE FIRST 
POINT). THE FFT GENERATES A SPECTRUM CONTAINING VALUES AT 
INCREMENTS OF OF IN FREQUENCY, WHERE DF = 1/(NX*DT). IF WE WERE 
TO TRANSFORM AN ARRAY OF NX COMPLEX DATA POINTS (WHOSE IMAGINARY 
VALUES WERE ZERO, SINCE WE ARE CONCERNED WITH A REAL SIGNAL) WE 
WOULD GENERATE NX COMPLEX SPECTRAL POINTS RANGING IN FREQUENCY 
FROM 2 TO (NX-1)*DF Hz. HOWEVER, THE POINTS FROM (NX/2+1)*DF Hz 
TO (NX-1)*DF Hz ARE MERELY THE COMPLEX CONJUGATES OF THE LOWER 
POINTS, AND CONTAIN NO ADDITIONAL INFORMATION. THIS IS BECAUSE 
THE NYQUIST RATE, OR THE HIGHEST UNALIASED FREQUENCY THAT CAN BE 
SAMPLED, IS 1/2*NX*DF. WHEN WE TRANSFORM NX REAL CLOSE-PACKED 
POINTS AND THEN EXTRACT THE REAL TRANSFORM, WE GET ONLY THE LOWER 
NX/2 COMPLEX SPECTRAL VALUES, WHICH ARE ALL THAT ARE NEEDED. 


ee we he te ee ee Ok 8 Oe OP OO we we 40 Oe Oe ee OF SE ee Oe ee ee we ee we we er oe ee en we ee SE OO we we te we we we ee we es we es we 
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5. NEXT, THE SPECTRUM IS CONVOLVED WITH A SHORT WINDOW FUNCTION. 
THE REAL AND IMAGINARY (EVEN AND ODD) POINTS ARE CONVOLVED 
SEPARATELY, SINCE THE WINDOW FUNCTION IS A SEQUENCE OF REAL 
NUMBERS. WINDOWING IS NECESSARY IF THE SIGNAL FILTERS WE ARE 
EFFECTIVELY SYNTHESIZING ARE TO HAVE A USEFUL SHAPE. WITHOUT 
WINDOWING WE WOULD FIND THAT EACH SPECTRAL FILTER WOULD HAVE A 
NARROW PASSBAND BUT SIGNIFICANT SIDELOBE RESPONSES. THAT IS, THE 
MAGNITUDE OF A SPECTRAL LINE X( 1) WOULD DEPEND NOT ONLY ON SIGNAL 
COMPONENTS NEAR 1*DF IN FREQUENCY, BUT ALSO SIGNIFICANTLY ON 
COMPONENTS AT OTHER FREQUENCIES AS WELL. WINDOWING BROADENS THE 
SHAPE OF THE PASSBAND OF EACH SPECTRAL FILTER BUT DECREASES THE 
SIDELOBE RESPONSES. THE AMOUNT OF BROADENING AND SUPPRESSION OF 
SIDELOBES DEPENDS ON THE WINDOW ORDER, OR THE LENGTH OF THE 
WINDOW FUNCTION WHICH IS CONVOLVED WITH THE RAW SPECTRUM. 


+ wee we we we te we we te ee ws es we we we 


6. FINALLY, THE MAGNITUDE OF THE SIGNAL AT EACH SPECTRAL 
FREQUENCY IS CALCULATED AS THE ROOT SUM OF THE SQUARES OF THE 
REAL AND IMAGINARY SPECTRAL COMPONENTS (IE, THE VECTOR MAGNITUDE 
OF EACH COMPLEX SPECTRAL POINT) AND SCALED TO MAKE] THE PROCESSING 
GAIN INDEPENDENT OF THE SIZE OF THE BLOCK LENGTH NX. WE COULD 
ALSO CALCULATE THE PHASE OF THE SPECTRAL POINTS BUT THIS 
INFORMATION ISN'T TOO VALUABLE FOR OUR USE. NOTE THAT THE 
HIGHEST MEANINGFUL SPECTRAL FREQUENCY DEPENDS ON THE ANTI- 
ALIASING LOW-PASS FILTER THAT YS USED WHEN THE SIGNAL IS SAMPLED. 
THERE IS NOT MUCH MEANING TO SPECTRAL POINTS ABOVE THE FILTER 
CUTOFF FREQUENCY, WHERE SIGNALS ARE ATTENUATED AND! FREQUENCY 
ALIASING BECOMES A PROBLEM. THUS, IN THIS PROGRAM WE DO NOT 
CALCULATE SPECTRAL COMPONENTS ABOVE ABOUT 38 Hz (THE NOMINAL 
CUTOFF FREQUENCY OF THE LOW-PASS FILTER}. 

3 

sTHIS PROGRAM IS WRITTEN IN BASIC FOR THE IBM PC COMPUTER. HOWEVER, IT 
31S EASILY ADAPTED TO OTHER MACHINES. FOUR THINGS TO WATCH OUT FOR WHEN 
s TRANSFERRING THE PROGRAM TO ANOTHER VERSION OF BASIC ARE: 


1. ALL ARRAYS START WITH INDEX @. THAT IS, THE ARRAY X(NX) RUNS 
FROM X(@) TO X(NX-1). ARRAY REFERENCES WILL NEED TO BE CHANGED 
IF THE PROGRAM IS TO BE USED ON MACHINES WHERE THE: FIRST ELEMENT 
OF ARRAYS HAS INDEX 1, 


we we me we ws we we we ne oe wn we we ce 


2. ALL VARIABLES WHOSE NAMES START WITH THE LETTERS "I" THROUGH 
"N" ARE IMPLICITLY INTEGERS (DEFINT I-N STATEMENT). THESE 
VARIABLES MAY NEED TO BE RENAMED. IT IS IMPORTANT) THAT LOOP 
COUNTERS AND ARRAY INDEX VARIABLES BE INTEGERS BECAUSE INTEGER 
ARITHMETIC (WHERE AVAILABLE) IS FASTER THAN FLOATING-POINT 
ARITHMETIC. 


3. THIS PROGRAM USES THE INTEGER DIVISION OPERATOR "\", THIS 
MAY BE CHANGED TO "/" IN OTHER VERSIONS OF BASIC. 


4, _ THIS PROGRAM USES LOGICAL OPERATORS (“NOT", "AND", “OR") ON 
INTEGER VARIABLES (NOT LOGICAL VARIABLES) IN THE BIT-REVERSAL 
ROUTINE. IF THESE OPERATORS ARE NOT AVAILABLE THE: ALTERNATE 
VERSION OF STATEMENTS 3988-3198 USING ONLY INTEGER ARITHMETIC 
MUST BE SUBSTITUTED. 
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3* * 
¢* CHOICE OF SAMPLING TIME AND BLOCK SIZE = 
3% * 


gE RARER ERIK RIE R HH RR HMR RWW RW ee 


3 

3THE SAMPLING TIME DT AND THE BLOCK SIZE NX TOGETHER DETERMINE THE 
;NUMBER OF SPECTRAL POINTS CALCULATED AND THEIR SPACING IN FREQUENCY. 
3THE CHOICE OF SAMPLING TIME ALSO AFFECTS THE FREQUENCIES AT WHICH POWER- 
sLINE HARMONIC INTERFERENCE THAT IS PASSED BY THE ANTI-ALIASING LOW-PASS 
sFILTER WILL APPEAR IN THE SPECTRUM. THE FOLLOWING DISCUSSION ASSUMES 
;THAT THE 38 Hz LOW-PASS FILTER WILL BE USED TO ATTENUATE HIGH FREQUENCY 
;SIGNAL COMPONENTS BEFORE SAMPLING. 


3 

;IF THE INPUT SIGNAL THAT IS SAMPLED CONTAINED ONLY FREQUENCIES BELOW 

338 Hz IT WOULD BE SUFFICIENT TO SAMPLE AT THE NYQUIST RATE OF 68 SAMPLES 
;PER SECOND. HOWEVER, SINCE THE FILTER DOES NOT HAVE INFINITE 
sATTENUATION ABOVE 38 Hz IT IS NECESSARY TO SAMPLE AT A SOMEWHAT HIGHER 
;RATE, AND THEN DISCARD THOSE SPECTRAL POINTS REPRESENTING SIGNALS ABOVE 
338 Hz. 


1. ONE CONSIDERATION IN CHOOSING THE SAMPLING RATE IS THE ACTUAL 
RATE AVAILABLE WITH A GIVEN A/D CONVERTER. WITH THE CMC BUSSTER 
D16 A/D CONVERTER, SAMPLING TIMES CAN BE AS SMALL AS &.8895 s. 
HOWEVER, IF MORE THAN ONE SIGNAL IS TO BE DIGITIZED AT ONE TIME, 
THE MINIMUM SAMPLING TIME IS 8.81 s PER CHANNEL, OR A SAMPLING 
RATE OF 188 SAMPLES/SECOND. THIS SEEMS LIKE A REASONABLE CHOICE. 


; 

3 

5 

3 

; 2. A SECOND CONSIDERATION IS THE EFFECT OF THE SAMPLING RATE ON 
: THE FREQUENCIES OF ALIASED POWER LINE HARMONICS. ONLY HARMONICS 
; AT 68 AND 12@ Hz ARE LIKELY TO BE A PROBLEM. WITH OT = @.@1 S, 

$ THE NYQUIST RATE IS 58 Hz, SO THE SPECTRUM WILL CONTAIN POINTS AT 
3 FREQUENCIES FROM @ TO 58 Hz. ANY 6# Hz SIGNAL THAT IS DIGITIZED 
; WILL APPEAR IN THE SPECTRUM AT 49 Hz, WHICH POINT WILL BE THROWN 
; OUT, SO 6@ Hz INTERFERENCE WON'T BE A PROBLEM. HOWEVER, 

; INTERFERENCE AT 128 Hz WILL APPEAR AT 28 Hz IN THE SPECTRUM, AND 
; THIS FACT MUST BE KEPT IN MIND WHEN ANALYZING THE DATA. (HIGHER 
3 HARMONICS APPEAR AS: 188 AT 28 Hz, 248 AT 48 Hz, 388 AT B Hz, 

H 368 AT 48 Hz, AND SO ON.) 

j 


3GIVEN A SAMPLING TIME DT = @.@1 s, WE CAN CALCULATE THE SPECTRAL SPACING 
sFOR DIFFERENT CHOICES OF BLOCK SIZE NX, AND WE FIND THE FOLLOWING: 


$ NX DF NUMBER OF POINTS <= 38 Hz 
3 - 64 1.5625 Hz 19 
3 128 @.7813 Hz 38 
; 256 @.3986 Hz 76 
$ 512 8.1953 Hz 153 
3 1924 @.8977 Hz 387 
$ 2048 2.8488 Hz 614 
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:* x 
3* WINDOWING * 
a * 
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§ 

sAFTER THE SPECTRUM OF THE INPUT SIGNAL HAS BEEN FOUND, IT !S CONVOLVED 
sWITH A SHORT WINDOW SEQUENCE TO IMPROVE THE SHAPE OF THE SYNTHESIZED 
sFILTERS. THOUGH THE SPECTRAL VALUES ARE COMPLEX, THE WINDOW SEQUENCE IS 
sA SEQUENCE OF REAL NUMBERS, SO THE CONVOLUTION IS DONE SEPARATELY FOR 
sTHE REAL AND IMAGINARY SEQUENCES IN THE SPECTRUM. THE CONVOLUTION IS OF 
sTHE FORM 

3 

$ XC J) <- W3*EX( 9-304XC 9439] +W2*0K( 9-2)4X( 9 42)) 

3 FWIF EXC gm LIFXC +193 tWO*XC]) 

3 

sWHERE X(j) IS THE REAL OR IMAGINARY SUBSEQUENCE OF THE SPECTRUM, AND THE 
sCONVOLUTION IS ILLUSTRATED FOR A 4TH ORDER WINDOW SEQUENCE. NOTE THAT 
sIF THE CONVOLUTION IS TO BE DONE IN-PLACE, A SMALL ARRAY ‘MUST BE USED TO 
sHOLD THE PREVIOUS VALUES OF X(3-3), X(y-2), X(g-1), and X(j9) FOR USE IN 
sSUBSEQUENT CONVOLUTIONS. 


3 : 

sTHE VALUES TO USE FOR DATA POINTS BEYOND THE ENDS OF THE SPECTRUM ARE 
sTHE COMPLEX CONJUGATES OF THE POINTS REFLECTED ABOUT THE ‘ENDS, SO THE 
sAUXILLIARY ARRAY OF PREVIOUS VALUES MUST BE INITIALIZED AS FOLLOWS: 

3 


; Re(-3) = Re(3) Im(-3) = -Im(3) 
; Re(-2) = Re(2) Im(-2) = -Im(2) 
3 Re{(-1) = Re(l) Im(-1) = -Im(1) 
3 


sPOINTS AFTER THE FAR END OF THE SPECTRUM ARE REFLECTED IN THE SAME WAY, 
+BUT SINCE WE WON'T CALCULATE SPECTRAL POINTS ABOVE 38 Hz THERE IS NO 
ae TO WINDOW POINTS THERE. 


THE FOLLOWING WINDOW SEQUENCES ARE THE MOST IMPORTANT FOR OUR USE. 
sTHESE ARE ALL MINIMUM-SIDELOBE WINDOWS WHICH SUPPRESS SIDELOBE RESPONSES 
1AS MUCH AS POSSIBLE FOR A GIVEN WINDOW ORDER, AND HAVE BEEN NORMALIZED 
sFOR UNITY GAIN. NOTE THAT THE ORDER-1 WINDOW aURE RE. CH TO NO WINDOW 

c 


;AT ALL. THE ORDER-2 WINDOW IS SIMILAR TO THE HAMMING FUNCTION. 

; ORDER 1 2 3 4 

; wa 1. 290g 1. e908 1.2980908 1.99980000 
; W1 g -8.42875 -8.5895613  -@.6727198 
; W2 g g @.8922278 | 9.1878524 
' W3 g g g -9.8146337 
$ 

i HIGHEST -13 dB -43 dB -72 dB -98 dB 

' SIDELOBE 

3 : 

; 3-dB WIDTH 8.89 DF 1.38 OF 1.61 DF 1.86 DF 

3 

; 6-dB WIDTH 1.20 OF 1.81 DF 2.25 DF 2.62 DF 


sFOR OUR USE, I THINK THE 3RD-ORDER WINDOW IS THE MOST APPROPRIATE, AND 
; THE PROGRAM USES THIS WINDOW FUNCTION. SINCE W3 = @ IN THIS CASE, THE 
sACTUAL CONVOLUTION CONTAINS ONLY 3 TERMS. SINCE W@ = 1, THE LAST TERM 
31S JUST X(3), WITHOUT A MULTIPLICATION. THUS THE WINDOWING OF EACH 
sPOINT INCLUDES 5 ADDITIONS AND 2 MULTIPLICATIONS FOR EACH OF THE REAL 
sAND IMAGINARY COMPONENTS. IF ONLY VERY SHORT BLOCKS (SAY NX = 64) ARE 
sTO BE PROCESSED, THE 2ND-ORDER WINDOW MIGHT BE MORE APPROPRIATE SINCE 
s;THE FILTERS ARE A BIT NARROWER IN THAT CASE. 

3 
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3" * 
ial PROGRAM * 
3” * 
gM OI ROO ROI IOI II TOTO TOR RI IORI IORI 
3 

;DEFINITIONS: 

18988 DEFINT I-N "VARIABLES STARTING “I"-"N" ARE INTEGERS 

1918 NX = 1824 "DATA BLOCK SIZE 

1928 NXM1 = NX-1 

1838 DIM X€NX) "INPUT DATA ARRAY X( 4), 1 = @, «22, NX-1 

1840 OT = 9.81 "INPUT SAMPLE TIME IN SECONDS 

1958 DF = 1./(DT*NX) "SPECTRAL LINE SPACING IN Hz 

1962 NMAX = 308./0DF "NUMBER OF HIGHEST SPECTRAL LINE <= 38 Hz 
1192 NF = NX\2 ‘FFT SIZE (NR OF COMPLEX POINTS) 

1118 MF = LOG(NF)/LOG(2.) "NF = 24MF 


1128 NFD2 = NFX\2 

1138 NFM1 = NF-1 

1148 DIM CA(NF) "COSINE AND SINE ARRAYS, 

1159 DIM SACNF)} " ANGLES FROM @ TO 188*(NX~-1)/NX 
116% PI2 = 6.2831853 


1229 Wl = -%.5859613 ‘WINDOW COEFFICIENTS 
1218 W2 = #.8922278 


A Wee WWW ee We We ee eR AI ee i II Cie a ie aie i we ee ie ae 


3 
sCALCULATE SINE AND COSINE ARRAYS: 


3 

sTHIS NEED BE DONE ONLY ONCE. NOTE THAT THE FFT ROUTINE USED SINES AND 
sCOSINES FROM & TO ¢18@ DEGREES IN STEPS OF 368/NF DEGREES, AND THE REAL 
3 TRANSFORM EXTRACTION ROUTINE USES SINES AND COSINES UP TO <9@ DEGREES IN 
s;STEPS OF 36@/NX DEGREES. THUS THE ODD NUMBERED POINTS IN THESE ARRAYS 
sFOR ANGLES ABOVE 98 DEGREES ARE NOT USED. IF MEMORY SPACE IS A PROBLEM, 
sTHE ARRAYS CAN BE CONDENSED AND ONLY POINTS FROM @ TO £18@ DEGREES IN 
3STEPS OF 36@/NF CALCULATED FOR USE BY THE FFT. IN THIS CASE, THE VALUES 
;USED BY THE REAL TRANSFORM EXTRACTION WILL HAVE TO BE CALCULATED AS 
sNEEDED. IT IS ALSO POSSIBLE TO USE A TABLE CONTAINING VALUES FOR ANGLES 
sONLY FROM @ TO 45 DEGREES, BUT THIS MAKES THE LOOKUP PROCEDURE SOMEWHAT 
sMORE COMPLICATED, INVOLVING A POSSIBLE CHANGE OF SIGN OR EXCHANGE OF SIN 
sAND COS VALUES DEPENDING ON THE SEMIQUADRANT OF THE ANGLE. 


1522 DA = PI2/NX "ANGLE INCREMENT 

1518 FOR J = # TO NFM1 "CALCULATE SINES AND COSINES 
1528 A = DA*) 

1538 CA(J) = COS(A) 

1548 SA( J) = SINCA} 


1558 NEXT 9 

JCC OI GIGIOIOIIOIIIGIDIIUICIOIOIOIOIOICIOIDISI IOUS IOI OIGIIDIGIUIOIIOI IOI OIGIIOI IO IO IDISIOOOIOIOI CIO 
3 

sSAMPLE INPUT DATA INTO ARRAY X: 

; THIS CODE MUST BE WRITTEN TO SAMPLE THE INPUT SIGNAL EVERY %.@1 SECOND 
sAND STUFF SUCCESSIVE VALUES INTO X. ONLY NX VALUES NEED BE TAKEN. 
sENTER HERE TO START PROCESSING A NEW BLOCK OF DATA. 


32828 CODE VET TO BE WRITTEN 
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3 
sBIT REVERSAL: 


$ 
3THE REAL COMPACT SAMPLES IN X ARE NOW TREATED AS A COMPLEX ARRAY OF NF 


3(REAL, IMAGINARY) PAIRS. 


sPRIOR TO THE FFT ROUTINE. 
3 
3THE FOLLOWING ROUTINE USES LOGICAL OPERATORS. 


PAIRS AT BIT-REVERSED ADDRESSES ARE NOW SWAPPED 


‘INITIAL BIT-REVERSED ADDRESS 
‘SCAN ALL ADDRESSES EXCEPT THE FIRST 
“INITIAL BIT MASK 


IF NOT(K AND L? GOTO 3978 'THIS BIT SET? 


' YES - CLEAR IT, 

* SHIFT MASK RIGHT, 

‘ AND CONTINUE 

‘BIT WAS CLEAR, SET IT 

‘SKIP IF ADDRESS SAME OR ALREADY DONE 
‘ARRAY ADDRESS OF LOWER REAL COMPONENT 
sARRAY ADDRESS OF UPPER REAL COMPONENT 
"SWAP REALS 


‘SWAP IMAGINARIES 


sTHE FOLLOWING VERSION USES ONLY INTEGER ARITHMETIC BUT IS A BIT LESS 


3929 K = 1 

3918 FOR J = 1 TO NFMI 
3828 L = NFD2 

3832 

3848 K = K AND NOT L 
3858 L = LN2 

3868 GOTO 3438 

3878 K = K ORL 

3888 IF K <= J GOTO 3198 
3898 JZ 2 Jtd 

3188 K2 = K+K 

3118 T = X(J2) 

3128 X(J2) = X(K2) 
31398 XCK2) = T 

3148 T = X(J2+1) 

3158 XCJ2+1) = X(K2+1) 
3189 XCK2+1) = T 

3199 NEXT J 

sEFFICIENT. 

3 

33928 FOR J = 1 TO NFMI 
33818 Ki = J 

33828 K = 2&8 

33838 FOR I = 1 TO MF 
33849 K2 = K1\2 

13858 K = (K-K2)*2+K1 
33969 Kl = K2 

33878 NEXT I 

33889 IF K <= J GOTO 3198 
33899 J2 = Jtd 

33199 K2 = K+K 

33118 T = X(J2) 

33128 XCJ2) = X(K2) 
33138 X(K2) = T 

33149 T= XCJ2+1) 

33158 XCJ2+1) = MCK241) 
33188 XCK2+1) = T 

33198 NEXT J 


“SCAN ALL ADDRESSES EXCEPT THE FIRST 
“INITIALIZE DIVISOR 

INITIALIZE BIT-REVERSED ADDRESS 
sBIT-REVERSE ALL MF BITS 


"SKIP IF ADDRESS SAME OR ALREADY DONE 
"ARRAY ADDRESS OF LOWER REAL COMPONENT 
sARRAY ADDRESS OF UPPER REAL COMPONENT 
"SWAP REALS 


"SWAP IMAGINARIES 
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3 
sFFT ROUTINE: 


;THIS ROUTINE USED DECIMATION IN TIME. THE COMPLEX DATA INPUT IS IN BIT- 
sREVERSED ORDER. 


1. THE OUTER LOOP MAKES MF PASSES THROUGH THE DATA. THE INITIAL 
NUMBER OF SUB-DFT'S IS 1, AND DOUBLES ON EACH PASS. THE INITIAL 
INCREMENT IN ANGLE BETWEEN SUB-DFT'S IS 188 DEGREES, AND HALVES 
ON EACH PASS. 


2. THE MIDOLE LOOP PROCESSES THE SUB-DFT'S. AT EACH SUB-OFT 
ANGLES INCREASE IN NATURAL ORDER AS @ (121); 8, 98 (I=2); @, 45, 
9@, 135 (I1=3); 8, 22.5, 45, ... (€I=4); AND SO ON. 


3. THE INNER LOOP CALCULATES THE BUTTERFLIES FOR A GIVEN ANGLE 
IN EACH SUB-DFT. EACH DECIMATION-IN-TIME BUTTERFLY HAS THE FORM 


Re(j) <- Rely d+ fRe(k )*C+#Im(k )*8] 
Im(3) <- Im€j)}+0ImCk )*C-Re(k >*S) 
Retk) <- Reg )-CRe(k *C+Im(k >*S] 
Im(k) <= Im€j)-CImCk )*C-Re(k)*S] 


WHERE C AND S ARE THE COSINE AND SINE OF THE ANGLE. 


Oe ae ee we we we te we we we we we we we we we we we a 


3288 LM = 1 ‘INITIAL NUMBER OF SUB-DFT'S = 1 

3219 LM2 = LM+LM 

3228 IDA = NF ‘INITIAL ANGLE INCREMENT = 189 DEG 
3239 FOR I = 1 TO MF ‘OUTER LOOP 

3246 IA = & ‘INITIAL SUB-OFT ANGLE = @ DEGREES 
3259 FOR L = 1 TO LM ‘MIDDLE LOOP 

3268 C = CACIA) "LOOK UP COS AND SIN FOR THIS ANGLE 
3278 S = SACTIA) 

3289 FOR J*L-1 TO NFM1 STEP LM2 ; INNER LOOP 

3298 JR = J+d "ARRAY ADDRESS OF LOWER REAL COMPONENT 
3388 JI = JR+1 "ARRAY ADDRESS OF LOWER IMAG COMPONENT 
3318 KR = JR+LM2 "ARRAY ADDRESS OF UPPER REAL COMPONENT 
3322 KI = KR+1 "ARRAY ADDRESS OF UPPER IMAG COMPONENT 
3338 TR @ XCKRI*C+XCKID*S ‘BUTTERFLY 

3348 TI = XCKID*C-XCKRINS 

3358 XCKR) = XCIRI-TR 

3369 XCKI) @ XCQI)-TI 

3378 XCIR) = XCIRD+TR 

3388 XCJT) = XCOT+TI 

3398 NEXT 9 "END OF INNER LOOP 

3498 IA = IA+IDA ‘BUMP ANGLE 

3419 NEXT L ‘END OF MIDDLE LOOP 

3349 LM = LM2 ‘DOUBLE NUMBER OF SUB-DFT'S 

3439 LM2 = LM2+LM2 

3442 IDA = IDAN\2 ‘HALVE ANGLE INCREMENT 

3458 NEXT I "BOTTOM OF OUTER LOOP 
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$ 
sEXTRACT CLOSE-PACKED REAL DATA TRANSFORM: 


$ 

s+NOW WE MUST EXTRACT THE SPECTRUM OF THE INPUT DATA AS IF IT HAD BEEN 
sCALCULATED FROM A SEQUENCE OF NX COMPLEX INPUT VALUES ALL OF WHOSE 

THE ALGORITHM FOR THIS IS AS FOLLOWS: 


s IMAGINARY COMPENENTS WERE ZERO. 


3 1. FOR POINTS 9 #= 1, 2, ..+, NF/2-1 CALCULATE 

' Rely) <- Ral y}+Relk )-[Re( 3 )-Rel(kd1*S+LImly t#Im(k °C 
3 Im(3) <- Im( gj }-ImCk }-CRel J )-Re(k I *C-CimGg d+Im(k)1*S 
: Ret{k) <- Rel y)+Relk )+CRe(y d-ReCk I*S-Cim(y d+Im¢k 1 *C 
; Im(k) <- -Im€3)+Im(k )-CRe ly )-Retk 1*C-CIm(y )+1m(k)1*S 
$ 

$ WHERE k = NF ~-3, A # g*368/NX DEGREES, C = COS(A), AND 

$ § = SINCA). 

$ 

3 WE ALSO HAVE THE SPECIAL POINTS: 

3 

3 Re(@) <- CRe€(H)+Im(9)] 

t Im(B) <- @ 

$ 

$ Re(NF/2) <~- Ret NF/2)¥*2 

3 Im(NF/2) <= -Im(NF/2)*2 

; AND (THOUGH THESE VALUES ARE NOT CALCULATED): 

; Re(NF) ¢- CRe(8)-Im(8)1*2 

3 Im(NF) IS UNKNOWN 


3599 X(H) = XCH)*2, 
3518 X(1) = 
3528 FOR J = 1 TO NFD2-1 


3532 JR = J+d 

3549 JI = UR+l 

3558 KR = NX-JR 

3562 KI = KR+tl 

3578 C = CA(J)- 

3589 S = SA(J) 

3598 RS = XCJRI+XCKR) 
3688 RD = XCJR)-XC( KR) 
3618 QS = XCJI)+X(KT) 
3628 QD = XC(JI)-X(KI> 
3638 TR = QS*C-RD*S 
3648 TI = RD*C+QS*S 
3658 XCJR) = RS+TR 
3669 X(J1T) = QS-TI 
3672 X(KR) = RS-TR 
3688 © X(KI) = -QS-TI 


3788 XCNF) = XONF)*2, 
3718 XCNF+1) = XCNF4#1)"2. 


"SPECIAL CASE FOR 


‘LOOP FOR GENERAL 
"ARRAY ADDRESS OF 
"ARRAY ADDRESS OF 
‘ARRAY ADDRESS OF 
‘ARRAY ADDRESS OF 


‘LOOK UP COS AND SIN FOR 


‘REAL SUM 
"REAL DIFF 
‘IMAG SUM 
‘IMAG DIFF 
‘TERM FOR REALS 
‘TERM FOR IMAGS 


"BOTTOM OF LOOP 
"SPECIAL CASE FOR 
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X(B) 


CASE 

LOWER REAL COMPONENT 
LOWER IMAG COMPONENT 
UPPER TNA COMPONENT 


UPPER IMAG COMPONENT 
'HIS ANGLE 
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; WINDOW SPECTRAL DATA TO IMPROVE FILTER SHAPE: 


;SEE THE DISCUSSION OF WINDOWING ABOVE. THIS ROUTINE IS FOR A 3RD-ORDER 
sWINDOW FUNCTION. 


3898 R2 2 X(4) ‘INITIALIZE PREVIOUS VALUE ARRAY, Re(2) 
3818 Rl = X(2) * Re(l) 

3828 RB = X(H) " Rel @) 

38398 Q2 = -X(5) '  Tm(2) 

3842 Ql = -X(3) '  Tm¢1) 

3859 Qg = -X(1) ‘ Imt(G) = B. 

3869 FOR J@ = @ TO NMAX*2 STEP 2 ‘WINDOW POINTS <= 38 Hz 

3878 Jl = JH+1 "ADDRESS OF IMAG COMPONENT 

3889 XC JD = CR24X(IMt4) MW24CRIT+KCIG+2)*W1tRD ‘REAL COMPONENT 
3898 XCIL) = CQ2+KCI1+4) *W24(Q1+XCI1+2))*Wl+Q9 ‘IMAG COMPONENT 
3998 R2 = Rl "SHIFT PREVIOUS ARRAY 

39198 Rl = R@ 

3928 RO = X(JH+2) ‘ AND INSERT NEW VALUE 

39398 a2 = Ql 

3949 Ql = Qs 

3958 Qe = X(J1+2) 

39698 NEXT J ‘END OF LOOP 


ge eK HR HR RRR KR RWW RK RK KKK RK 


3 
;CALCULATE MAGNITUDES: 


; THE MAGNITUDE OF EACH SPECTRAL VECTOR IS CALCULATED AS 
$ 
3 MAGNITUDE( 3) = SQRT[Re( 3 )*Rel 3 )+Im(y>*Im( 9] *SF 


5 
sWHERE SF = 1/NX IS A SCALING FACTOR, USED FOR 2 REASONS: 


1. THE SCALING FACTOR SETS THE PROCESSING GAIN. WITH SF = 1/NX, 
THE GAIN FOR A SINGLE-COMPONENT SIGNAL IS 1. THAT IS, IF THE 
INPUT SIGNAL IS A SINE-WAVE AT FREQUENCY n*DF WITH PEAK AMPLITUDE 
1.8, THE SPECTRUM AT FREQUENCY n*DF WILL ALSO BE 1.8. 


2. THE SCALING FACTOR MAKES THE PROCESSING GAIN INDEPENDENT OF 
THE BLOCK SIZE NX, SO THE SAME SIGNAL CAN BE PROCESSED WITH 
DIFFERENT BLOCK SIZES IN ORDER TO GENERATE DIFFERENT SPECTRA WITH 
VARYING FREQUENCY RESOLUTION, BUT THE MAGNITUDE OF A NARROW-BAND 
SIGNAL COMPONENT WILL BE THE SAME IN EACH SPECTRUM. 


we es ee es ws se we we es we we 


3 

sNOTE THAT MAGNITUDES ARE ONLY CALCULATED FOR SPECTRAL POINTS <= 38 Hz. 
sMAGNITUDES ARE STORED AT THE BOTTOM OF THE ARRAY X, IN THE FIRST NMAX 
;POINTS. 


3 

;THIS ROUTINE CAN BE INCLUDED IN THE WINDOWING ROUTINE ABOVE WITH A SMALL 
s INCREASE IN PROGRAM EFFICIENCY, BUT IT IS SHOWN SEPARATELY HERE FOR 
sCLARITY. 


4929 SF = 1./NX "SCALE FACTOR 

491g FOR J = # TO NMAX "PROCESS POINTS <= 398 Hz 

4828 J2 = JtJd "ADDRESS OF REAL COMPONENT 

4938 K2 = J2+1 "ADDRESS OF IMAG COMPONENT 

4849 XCF) = SQRTCXCI2)*XC9Z9+X(K29"*X(K2)9*SF ‘CALCULATE MAGNITUDE 
49SD NEXT J ‘BOTTOM OF LOOP 


;END OF PROGRAM. 
3 
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3* * 
Bal TESTING THE PROGRAM * 
7% * 


He Wee We We We We ee We We Wee Wee ei eee ee ee a 


;THE BEST WAY TO TEST THE PROGRAM IS TO GENERATE A SYNTHETIC INPUT ARRAY 
sOF DATA POINTS CONTAINING SIGNALS AT KNOWN FREQUENCIES, AND PROCESS THEM 
3AND CHECK THAT THE RESULTING SPECTRUM IS CORRECT. TO GENERATE A SIGNAL 
; ARRAY WHICH CONTAINS A SIGNAL AT A FREQUENCY Q*DF USE TH§ FOLLOWING CODE 
sBEFORE THE BIT-REVERSAL ROUTINE: 


2288 Q*... ‘FREQUENCY IN UNITS OF DF 
2818 OP = PI2*Q/NX "PHASE INCREMENT PER SAMPLE 
2828 FOR I = @ TO NXM1 "DATA GENERATION LOOP 

2838 XCT) = SINCOP*T ) “CALCULATE A DATA POINT 
20848 NEXT I "END OF LOOP 


3Q IS A REAL NUMBER, AND SHOULD BE CHOSEN SO THAT @ <= Q <= NX/2. NOTE 
sTHAT THE ARRAY X WILL CONTAIN @ CYCLES OF A SINEWAVE SIGNAL. 


sAFTER PROCESSING, THE SPECTRUM SHOULD CONTAIN ONLY COMPONENTS NEAR THE 
+FREQUENCY Q*DF. THAT IS, SPECTRAL LINES X(3) SHOULD BE SMALL EXCEPT FOR 
3J NEAR Q@. IF Q@ IS AN INTEGER, THE SPECTRAL LINE X(Q) SHOULD BE 1 AND 
sALL THE OTHERS SHOULD BE @ (OR OF ORDER 1E-6, DEPENDING ON THE ROUND-OFF 
sERRORS IN THE PROCESSING). USING NON~INTEGER VALUES OF Q WILL RESULT IN 
3A SPECTRUM WITH A CLUSTER OF RESPONSES NEAR X{INT(Q)) BUT LINES FAR AWAY 
sFROM Q SHOULD BE SMALL. THE FALLOFF IN RESPONSE IN THE SPECTRUM AS ONE 
sMOVES AWAY FROM THE LINE XCINT(Q)) DEPENDS ON THE EFFECTIVENESS OF THE 
sWINDOW FUNCTION, AND THE SHAPE OF THE SYNTHESIZED FILTER RESPONSE CAN BE 
sPLOTTED BY USING VARIOUS VALUES OF Q. 


5 

sNOTE THAT IF NMAX ¢ Q@ <= NX/2, THE SIGNAL GENERATED LIES IN THE UPPER 
sPART OF THE SPECTRUM, WHICH IS THROWN AWAY (DON'T EXAMINE SPECTRAL LINES 
3X€j) WHERE gy > NMAX, AS THEY HAVEN'T BEEN CALCULATED). ALSO, IF 

3NX/2 ¢ Q ¢ NX, THE SIGNAL WILL BE ALIASED IN THE SPECTRUM, APPEARING AT 
3A FREQUENCY OF NX/2-Q. 


; é 
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ULF TREE POTENTIALS AND GEOMAGNETIC PULSATIONS (U) 
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HIGH-SENSITIVITY measurements of  ultra-low-frequency 
(ULF; frequencies less than 5 Hz) geomagnetic pulsations’ 
usually require elaborate receiving antennas ranging from 
large air-cored coils’? through multi-turn steel, or mumetal- 
cored solenoids*"* to small superconducting loops immersed 
in liquid helium". Pairs of electrodes inserted in the ground 
have also been used as antennas’*. The need for a large 
spacing between the electrodes (varying from hundreds to 
thousands of metres) and the difficulty of calibrating the 
measurements absolutely have resulted in the almost univer- 
sal use of the more compact and easily calibrated coil-type 
antennas in recent years. | describe here a new method for 
measuring ULF geomagnetic pulsations, which requires a 
minimum of elaborate equipment. The method is based on 
the use of trees, or, more specifically, on the use of pairs of 
electrodes inserted into trees, as ULF receiving antennas. 

There are several reasons that this new method of 
measurement may be of interest. The equipment is simple 
and thus the method could lead to more widespread observa- 
tions of ULF geomagnetic pulsation phenomena. The 
method of measurement also provides new information 
about tree potentials, that is, it shows that some, and 
perhaps all, of the ULF components of these potentials are 
induced by ULF geomagnetic field fluctuations and do not 
originate in the trees themselves. Finally, although it is not 
clear at present what effect induced ULF electric fields may 
have on the growth and other vital processes in a tree, the 
link between these ULF electric fields and geomagnetic 
field fluctuations suggests that some environment-related 
changes in trees could also be influenced by changes in 
geomagnetic activity. These changes may have a natural 
origin (for example, the changes that occur during a solar 
cycle*) or they may be caused by a variety of human activi- 
ties (by modern d.c.-powered mass transit systems, which 
can produce large amplitude ULF electromagnetic fields'’). 

The ULF measurements reported here were stimulated 
by the work of Burr on relatively steady-state tree poten- 
tials''’. Burr recorded these potentials for more than a 
decade using a pair of specially-designed non-polarisable 
electrodes inserted in the cambium of an unspecified tree 
(which was probably a maple). The electrodes were about a 
metre apart along the long axis of the tree and Burr 
observed diurnal, 27-d, and seasonal variations, as well as a 
suggestion of a correlation with sunspot activity, in their 
potential difference. 

Most of Burr's observations were at frequencies far 
below the frequency range for ULF geomagnetic pulsations. 
One series of measurement obtained, however, during an 
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ULF TREE POTENTIALS AND GEOMAGNETIC PULSATIONS 


electrical storm suggested that ULF variations of tree 
potentials might occur on occasion. I therefore began a 
search for. variations with frequencies predominantly in the 
Pcl geomagnetic pulsation range (0.2-5Hz). These 
frequencies correspond approximately to the delta regime 
for human brain waves. 

The measurements were made using a large native oak, 
Quercus lobata, that was located near conventional ULF 
recording equipment at a site on the Stanford University 
campus. This latter equipment uses 20,000 turn steel-cored 
solenoids as ULF antennas and it operated continuously 
throughout the interval during which the tree measurements 
were made. Thus, simultaneous measurements of ULF geo- 
magnetic pulsations using both conventional loop antennas 
and a tree ‘antenna’ were obtained at the one location. _ 

Two steel nails were used as electrodes. Following Burr’s 
configuration, they. were inserted about 0.05 m into the tree 
along the long axis, with a spacing of 0.76m. The lower 
electrode was approximately 1 m above the ground, and the 
two electrodes faced toward the geomagnetic west. Because 
the tree was not completely vertical, a line joining the two 
electrodes would have been inclined approximately 20° 
toward the geomagnetic east. The diameter of the tree mid- 
way between the two electrodes was 0.65 m. 

A resistance of about 5 k& was typically observed between 
the electrodes, increasing to about 10 kQ if polarisation was 
allowed to occur. A d.c. potential difference was also 
observed that varied from day to day but whose absolute 
value was usually in the range 10 to 100 mV, with the upper 
electrode positive. The electrodes were connected to a low- 
frequency high-gain amplifier through an RC filter (R= 
22 M2, C=50 uF). The amplifier was usually set for 50 db 
gain, and its output was filtered (0.02-7 Hz) before being 
recorded, generally without additional amplification, on a 
chart record and on analog magnetic tape. 

The ULF signals measured by this system were 
undoubtedly induced in the tree ‘antenna’ and not in the 
shielded cabling between the electrodes and the recording 
system: when the electrodes were disconnected from the 
tree and connected to an equivalent 5k resistor, without 
any other change in the wiring or configuration of the 
system, only a steady low level of white noise (typical resistor 
thermal noise) was observed. 

Similarities between the ULF signals recorded conven- 
tionally and with the tree ‘antenna’ were immediately 
apparent on the chart records. More detailed analysis con- 
firmed that Pcl pulsation events recorded by the two 
systems were very nearly identical in all their important 
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characteristics. Figure !, for example, shows spectrograms 
of a sequence of four Pc! pulsation events that occurred 
during the interval 1200 to 1500 ur on 17 January 1976, 
and which were received by the tree ‘antenna’ (a) and the 
conventional north-south solenoid antenna (b). With the 
exception of a lower signal-to-noise ratio for the tree 
measurements, the two Pc1 pulsation records are closely 
alike. It will also be noticed that the lower frequency Pc 2/ 
Pc 3 geomagnetic activity (frequencies in the range 0.02 to 
0.2 Hz) is recorded similarly by both systems. The amplitude 
of the ULF pulsations in the tree potentials is very small. 
For the Pc! pulsations shown in Fig. 1, the maximum 
amplitude of the potential fluctuations was about 0.1 mV. 
The nearly identical occurrence and spectral character- 
istics of ULF events measured by the tree electrodes and 
by the conventional ULF equipment indicated that the tree 
potentials were largely induced by ULF time variations of 
the geomagnetic field. To investigate this possibility, a 
portable planar search: coil powered by a 1! Hz signal 
generator was moved around the tree near the electrodes. 
It was found that a 1 Hz oscillation of the potential dif- 
ference between the tree electrodes was produced only when 
the search coil was orientated with its moment vector in 
the north-south direction. When the two electrodes were 
moved to the north face of the tree, a response from the 
electrodes could be obtained only when the search coil 
moment vector was orientated in the east-west direction. 
These results, and the observations of natural Pci pulsa- 
tions, can possibly be best understood by considering. the 
tree/electrode pair combination to form a collection of 
conducting loop antennas in which e.m.fs may be induced 
by magnetic field fluctuations in the appropriate direction. 
The conducting paths are provided by the conducting 
material of the tree (and the cambium in particular'’), and, 
for field fluctuations in a particular direction, the area of 


2.0-.7 


AEE OS 


Frequency (Hz) 


the relevant loop antenna is defined by the intersection of 
the tree with a vertical plane perpendicular to the particular 
field direction and passing through the two electrodes. Thus, 
in the measurements reported here, the Pc | pulsation events 
observed in the tree potentials were produced by Pc 1 pulsa- 
tions of the north-south component of the geomagnetic 
field. 

Further tests showed that the tree potentials could only 
be detected in a living tree. Thus, when a tree dies, the 
potentials gradually disappear as the wood dries and loses 
its conductivity. 

In conclusion, measurements with tree electrodes show 
that trees may be used as ‘antennas’ to detect ULF geo- 
magnetic pulsations. The measurements also show that ULF 
tree potentials are largely produced by ULF fluctuations of 
the geomagnetic field (the remaining component of the 
potentials is probably thermal noise). Presman’” noted that 
electromagnetic fields usually have an adverse effect on 
living processes. If the ULF geomagnetic pulsations have 
any adverse effect on the growth of trees (and, as we have 
seen, they must induce electric currents in the living 
material) these effects could possibly be observed in tree 
ting data. Pel geomagnetic pulsation occurrences vary 
markedly over a solar cycle’ and thus, if these particular 
pulsations effect tree growth, a solar cycle in tree ring data 
could occur. LaMarche and Fritts’ searched unsuccessfully 
for a relation between tree ring data and sunspot numbers. 
The phase of the Pc I pulsation solar cycle, however, differs 
by several years from the sunspot cycle and, assuming the 
two cycles affect tree ring data, they may tend to obscure 
each other's effects. Furthermore, other geomagnetic pulsa- 
tions and higher-frequency electromagnetic signals have 
their own cycles of occurrence, and their effects on tree 
ring formation, if any, could add further to the complexity 
of the tree ring data, Studies of these possible effects are 


Fig. 1 Spectrograms of a_ series 
of Pc 1 geomagnetic pulsation 
events recorded at Stanford, Cali- 
fornia, using tree potentials (a) 
and a conventional solenoid an- 
tenna (5). Short intervals of a 
1 Hz calibration signal appear 
at the start of each hour. The 
vertical lines in the upper spectro- 
gram are caused either by local 
electromagnetic transients or by 
natural sferics; similar lines occur 
in the lower spectrogram, but 
they are not as obvious because 
the background noise is com- 
paratively suppressed. 
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desirable, because the tree ring data could provide a unique 
record of past ULF and higher-frequency geomagnetic 
activity. 
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FORMAT FOR SATELLITE BROADCAST OF SPACE ENVIRONMENT SERVICES (U) 
(15 February 1984) 


There are seven sixty character Iines. Each is followed by a carriage 
return and a line feed. 
1234567890 1234567890 1234567890 1 234567890 1 2345678901 234567890 


In the examples of lines, CAPITAL letters are caption Information 
that appears as it does in the broadcast. Small letters represent 
data that are explained below. The : Indicates the 61st space 

(: Is not broadcast). Any data not avallable for broadcast will be 
replaced by ***, 


LINE ONE = GOES X=RAY AND TOTAL MAGNETIC FIELD 
(1 min values broadcast each minute) 


doy hhmm GOES6 Ix.x GOES5 Ix.x BOU TF ggggg GAMMAS 
347 2153 GOES6 B2.0 GOES5 B2.0 BOU TF 55397 GAMMAS 


where doy = day of year, such as 347 

hhmm = hours and minutes, such as 2153 

Ix.x = xeray level,(1 min avg), expressed with the B,C,M,X 
scale, where B = 1077 watts / meter2, C = 1076 
M= 10-5 and X = 10 -4, 
X-rays are measured in a 1-8 angstrom channel on the 
GOES geosynchronous satellites. An Increase In X-ray 
levels Indicates solar flare activity. 
In the example B2.0 = 2.0 * 107-7 watts / meter2. 


ggggg = Boulder mag fleld in gammas,(1 min avg), 
such as 55397 
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LINE TWO....GOES GEOSYNCHRONOUS MAGNETIC DATA 
(1 minute values broadcast each minute) 


GOES5 HPsgggg HEsgggg HNsgggg GOES6 HPsgggg HEsgggg HNsgggg : 
GOES5 HP+0123 HE-0023 HN+0033 GOES6 HP+0123 HE-0023 HN+0033 


where sgggg = sign (+ for positive, - for negative) and value 
(In gammas) of the magnetic fleld components at 
geosynchronous altitudes (approximately 22,000 
miles). These are measured by the GOES 5 (W075) 
and GOES 6 (W135) satellites. These components are 
HP (parallel! to the Earth's rotation axis,+ Is 
north), HE (directed earthward, + is down) and HN 
(perpendicular to the other two, le east or west, 
+ Is west). Example +0123. 


LINE THREE.... INDICES 
Minute 0 


IMEV p.pE+p 10MEV p.pE+p 100MEV p.pE+p PCA aa.aDB NTRN nnnn : 
1MEV 3.1E+O 10MEV 2.5E+0 100MEV 1.2E+O PCA 00.0DB NTRN 0451 


number of protons / (sec*cm2*ster) 

at greater than 1MEV, 10MEV,or 100MEV energies 

(5 min avg), such as 3.1E+0 ; 

DB of polar cap absorption at 30 MHZ as 

measured at Thule, (15 Min value), such as 00.0 
nnnn = Thule neutron monitor data In millivolts (15 min 
value). The background level is 460 to 470 mv. 
Significant deviations are + or - 10% or less than 
410 mv and greater than 510 mv. A sharp Increase 
that lasts for a few tens of minutes to hours is a 
Ground Leve! Event that Is caused by relativistic 
particles associated with solar flares. A sig- 
nificant decrease that lasts for several hours 

to several days Is a Forbush Decrease. This 
decrease Is the scattering of background cosmic 
rays due to plasma density enhancements and 
associated magnetic field enhancements In the near 
earth solar wind. These enhancements may or may 
not be solar flare related. 


where p.pEp 


aa.a 
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doy hh-hh K k/ggg BO Knnnn nonn PL Knnnn nnnn 10CM rrr 
347 18-21 K 3/033 BO K3334 333* PL K4334 433* 10CM 104 


where doy day of year (current day), such as 347 


hh=hh beginning and ending hours (in Universal 
Time or Z) for the 3 hr period of the K 
index to follow, such as 18=21 
k/ggg = Boulder K value / Gammas of deflection from the 


quiet day curve, such as 3/033 
mann nnnn = K Indices determined so far for this day (Boulder 

or planetary). The first value Is for 002-032, 
and each Index is for the 3 Hr perlod fol lowing 
the K index preceding It. Missing values wil] 
have an *, Example 3324 333* or 4334 433% 

rer = 10 CM solar radio flux In Solar Flux Units, 
(10°22 watts / (meter2*Hz)) as measured by 
Ottawa at approximately 1700Z each day. This value Is 
for the doy indicated at the start of the line. 
Until this value Is available *** will be broadcast. 
Example 104 


Minute 2 


doy BO Aaaa Knnnn nnnn PL Aaaa Knnnn nnnn TEDttt 
346 BO A017 K3334 4323 PL A024 K4333 4434 TED*** 


where doy = day of year (yesterday), such as 346 


aaa A Index (Boulder or Planetary) for yesterday, 
such as 017, or 024 
nnnn nann = K indices (Boulder or Planetary) for yesterday, 
(starting at 00-03Z) missing values will contain *. 
Example 3334 4323 or 4333 4434 
ttt = Total Energy Deposition, an estimate of the total 


power being dissipated In one hemisphere's auroral 
atmosphere by precipitating particles. The units 
are gigawatts. This Is a dally average for the doy 
Indicated. Estimates are made using particle data 
from Tiros satellites. Dally averages vary from 
near 005 on a quiet day up to about 100 on a very 
active day. If this data Is unavailable *** will 
be used. Example *** 
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Minute 3 


doy 10CM rrr 90 DAY MEAN mmm SSN sss XRAY BOD !x.x 
346 10CM 104 90 DAY MEAN 110 SSN 090 XRAY BGD B2.5 


where doy = day of year (yesterday), Such as 346 


ere 10 CM solar radio flux In Solar Flux Units, 
for 1700Z on doy Indicated, such as 104 

mmm = 90 day mean of the 10 CM Solar radio flux 
through yesterday, such as 110 

sss = Sun Spot Number for yesterday, such as 090 

Ix.x = Xray Background level for yesterday expressed on 
the B,C,M,X scale, where B = 1077 watts / meter2 
C= 10-6, M= 1075, and X = 1074.Iin the example, 
B2.5 = 2.5 * 107-7 watts / meter2. 

Minute 4 
doy GOES5 (W075) FLUENCE IMEV f.fE+f 10 MEV f.fEtf : 


346 GOES5 (WO75) FLUENCE 1MEV 2.8E+4 10 MEV 5.2E+3 


where doy = day of year (yesterday), such as 346 
f.fE+f = dally total fluence (protons / (cm@*day*ster)) 
greater than 1MEV or 10MEV (as indicated) expressed 
as an exponential number, such as 2.8E+4, or 5.2E+3 


LINE FOUR....HOURLY !ONOSPHERIC DATA 
Each 24 hour set of data will take 3 minutes to transmit 
(8 hours of data per |Iine). EMPTY Indicates that no data 
is being transmitted. 


sss ddddd doy hh nn.n nn.n nn.n van.n on.n van.n Aner Ann H 
BOU TEC 347 16 09.5 14.0 20.3 22.8 22.3 22.8 *#HE HEX 
where sss = station Identifier such as BOU: 
ddddd = data Identifier, such as TEC (Total Electron 


Count) ,FOF2 (Critical Frequency of the F2 Layer), 
M3000 (Optimum Frequency for a Single Hop 3000 KM 
Transmission), FMIN (Minimum Usable Frequency), 
or FOES (Critical Frequency of the Sporadic 

E Layer). 

hour of first data group, such as: 16 

data group. There are elght consecutive data 
groups. Missing data, or data not yet 

determined will have ****, Data for TEC 

is In units of 10+16 electrons / meter2, 

for the others It Is MHZ. Example 09.5 


hh 
nn.n 
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LINE FIVE....FORECASTS AND WARNINGS 
Minute 0 


doy-doy Mnn/nn/nn Xnn/nn/nn Pan/nn/nn 10CMerr/ere/err/reer 
348-350 MO1/01/01 X01/01/01 P01/01/01 10CM102/100/098/*** 


where doy-doy = day of year, first and last days of the forecast 
nn/nn/nn = percent probability for occurrence of activity 
on the first, second, and last day of the forecast. 
These apply to M flares, X flares, or Proton events 
as indicated. Example 01/01/01 
rrr/rrr/rer/rree = forecast 10 CM solar radio flux for the 3 forecast 
days, such as 102/100/098/*** The last group of 


rer is for the last doy plus one and will be *** 
until approx 18z on the first doy. At this time 
values will be added. 


Minute 1 


doy-doy AFaaa/aaa/aaa APaaa/aaa/aaa Kkkkk kkkk kkkk kkkk 
348-350 AFO18/015/015 AP020/020/020 K**** **33 3444 3333 


where doy-doy = day of year, first and last days of the forecast 
aaa/aaa/aaa = forecast A values for the forecast pertod,either 
AF, or AP as Indicated. Example 018/015/015, or 


20/20/20 
kkkk kkkk kkkk = Forecast KP values (without the 0,+,or-). 
kkkk The first of these 16 consecutive K values 


is for 00Z through 03Z on the first doy IIsted. 
Once a forecast K value has been measured It will 
appear in Iine two and the forecast value In this 


line will be replaced by an *, Example ¥*** ###3 
3444 3334. At approximately 18Z the remaining 
values will be updated with a new forecast. 
Minute 2 
PVITVITGITE... 
MAGALERT MINOR 
_where IIItl...= Description of forecast magnetic activity 
This will usually be: which Indicates: 
MAGQUIET only sporadic weak 


geomagnetic activity 
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MAGALERT JJ/KK or 
MAGALERT MINOR JJ/KK 


Increase in or continued 
high magnetic activity 
(JJ/KK are the first and 
fast calendar days of 


activity) 
MAGNIL end of active perlod or beginning 
of period of very low activity 
(used to end alerts) 
Minute 3 
PEEUTTIUETT ee 
PROTONNIL 
where II! 11...= Description of forecast proton activity 


This will usually be: 


PROTONQUIET 


which tndicates: 


no proton enhancements measured 
In, or forecast for the near 
earth environment. 


PROTON FLARE ALERT JJ/KK QXXYY 


Protons expected near earth 
because the solar region 
located In quadrant Q (1=NE, 
2=SE, 3=SW, 4=NW), Central 
meridian distance XX, latitude 
YY, has a 15% or greater 
probabilty of a proton 

flare from days JJ to KK 


PROTON ARRIVAL ALERT PP/FF HHmm 


PROTONALERT JJ/KK 


PROTONNIL 


Protons expected near earth 
on day PP from flare which 
occurred on day FF at HHmm 


Protons expected to meet the 
SESC alert criteria (see line 6 
ALERTS) In the near earth 
environment from days JJ to KK 


end of active period or begining 
of period of very low activity 
(used to end alerts) 
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LINE FIVE....FORECASTS AND WARNINGS (continued) 


Minute 4 
PUGTtdtdtbid... : 
SOLQUIET 
where [IHI11...= Description of forecast solar activity 
This will usually be: which Indicates 
SOLQUIET only sporadic weak solar activity, 


(class C or smaller flares) 


SOLALERT JJ/KK M Class flares, or an 
Increase in M level 
activity from days JJ 
to KK 


MAJOR FLARE ALERT JJ/KK QXXYY 
region located at QXXYY 
(see minute 3) has a 15% 
or greater prob of X class 
flares from days JJ to KK 


SOLNIL end of active period or begining 
of period of very low activity 
(used to end alerts) 


LINE SIX....ALERTS 
(rotates through all activated alerts) 
PUtdEdtdd. : 
ALERT NIL 
where IIIII...= descriptive text of alert conditions that are met 
These are The current alerts 
O.ALERT NIL (when no alerts are turned on) 
1.245 MHZ BURST > 100 SFU 
2.245 MHZ RADIO NOISE STORM > 5 TIMES BKGND 
3.10 CM FLARE > 100% ABOVE BKGND 
4.TYPE !! RADIO BURST 
5.TYPE IV RADIO BURST 
6.SUSPECTED PROTON FLARE 


7.PROTON EVENT 10 PART/CM2*SEC*STER > 10 MEV 
8 PROTON EVENT 100 PART/CM2*SEC*STER > 100 MEV 
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(the numbers are 9.XRAY >= M5 

transmitted with 10.XRAY >= X1 

the alerts) 11.A INDEX >= 20 (BOU) (forecast or observed) 
12.A INDEX >= 30 (BOU) " " " 
13.A INDEX >= 50 (BOU) 7 " " 
14.K INDEX >= 4 (BOU) (fest, obsvd,or in progress) 
15.K INDEX >= 5 (BOU) " te ta 
16.K INDEX >= 6 (BOU) " eo aS ee i 
17.K INDEX >= 5-5 (BOU) " Ae " 
18.SUDDEN STORM COMMENCEMENT (forecast or observed) 
19. STRATWARM 
20.SST ALERT: RADIATION >= 10 MILLIREMS (forecast or 

observed) 


LINE SEVEN...MESSAGES 
PITUUEITI. i 


where II111...= unformatted message with up to 60 characters wil 
sometimes be broadcast In thls line. Messages that 
will take more than one line will be broadcast 
sequentially over several minutes. These sequential 
messages will start with MMMM at the beginning of 
the first line. NNNN at the beginning of a IIne 
will Indicate that the message Is complete. 


In the future routine products will be transmitted 
ustng this line at preselected times. A tentative 
schedule of these products fol lows: 


PRODUCT , TIME 
Current message schedule 00Z 
SMM (Solar Maximum Mission) 01Z,07Z,13Z,19Z 


Observing Pian 

Region Report (body of 0030Z report) 02Z,082,14Z,20Z 
27 Day Forecast 03Z,09Z,15Z,21Z 
SGAS 

(Solar Geophysical Activity Summary) 04Z,10Z,16Z,222 
(body of 0245Z report) 

SDF (Report of Solar and Geophysical 
Activity..formerly Space Disturbance 052Z,11Z,17Z,23Z 


Forecast) (text portion of 22002: 
report) 


GEOALERT (body of 0330Z GEOALERT) 06Z,12Z,18Z 
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